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The formation of the beam halo in charged particle accelerators is studied in the framework of a 
stochastic-hydrodynamic model for the collective motion of the particle beam. In such a stochastic- 
hydrodynamic theory the density and the phase of the charged beam obey a set of coupled nonlinear 
hydrodynamic equations with explicit time-reversal invariance. This leads to a linearized theory 
that describes the collective dynamics of the beam in terms of a classical Schrodinger equation. 
Taking into account space-charge effects, we derive a set of coupled nonlinear hydrodynamic equa- 
tions. These equations define a collective dynamics of self-interacting systems much in the same 
spirit as in the Gross-Pitaevskii and Landau-Ginzburg theories of the collective dynamics for in- 
teracting quantum many-body systems. Self-consistent solutions of the dynamical equations lead 
to quasi-stationary beam configurations with enhanced transverse dispersion and transverse emit- 
tance growth. In the limit of a frozen space-charge core it is then possible to determine and study 
the properties of stationary, stable core-plus-halo beam distributions. In this scheme the possible 
reproduction of the halo after its elimination is a consequence of the stationarity of the transverse 
distribution which plays the role of an attractor for every other distribution. 

PACS numbers: 02.50.Ey, 05.40.-a, 29.27.Bd, 41.85.Ew. 



I. INTRODUCTION 



High intensity beams of charged particles, in par- 
ticular in linacs, have been proposed in recent years 
for a wide variety of accelerator-related applications: 
drivers for sources of neutron spallation; production of 
tritium; transmutation of radioactive wastes to species 
with shorter lifetimes; heavy ion drivers for fusion-based 
production of thermonuclear energy; production of ra- 
dioactive isotopes for medical use, and so on. In all these 
cases it is very important to keep at a low level the beam 
losses to the wall of the beam pipe, since even small frac- 
tional losses in a high current machine can cause exceed- 
ingly high levels of radioactivation. It is now widely be- 
lieved that one of the relevant mechanisms for these losses 
is the formation of a low intensity halo relatively far from 
the core of the beam. These halos have been either di- 
rectly observed Q or inferred from experiments Q, and 
have also been predicted from extensive numerical simu- 
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lations HSHESEHIHITIl. 0ften in these stud " 

ies the particle-core model has been used to understand 
the halo formation and its extent. That notwithstand- 
ing, it is widely believed that for the next generation of 
high intensity machines it is necessary to obtain a more 
quantitative understanding not only of the physics of the 
halo, but also of the beam transverse distribution in gen- 
eral [IJilIjJllJ]- In fact "because there is not a consensus 
about its definition, halo remains an imprecise term" [Tfij 
so that several proposals have been put forward for its 
description. 

Charged particle beams are usually described in terms 
of classical deterministic dynamical systems. The stan- 
dard model is that of a collisionless plasma where the 
corresponding dynamics is embodied in a suitable phase 
space (see for example 0). Generally speaking, the 
beams are described in a comoving reference system so 
that we can confine ourselves to a non rclativistic set- 
ting. In this framework, the formation of the halo has 
been studied mainly by means of the particle-core model 
initially proposed by Gluckstern, and the simulations 
show that the dynamical instabilities due to a paramet- 
ric resonance can allow ions to escape from the core 

in ma Hi nana. 

In the present paper we propose and develop a different 
approach. We introduce a model for the formation of the 
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halo based on the idea that the trajectories of the charged 
particles are sample paths of a stochastic process, rather 
than the usual deterministic (differentiable) trajectories. 
A random collisional plasma described in terms of con- 
figurational stochastic processes seems a realistic model 
for the dynamics taking place in particle accelerators. 
Moreover, it can explain the rare escape of particles from 
a quasi-stable beam core, since it takes into account, sta- 
tistically, the complicated inter-particle interactions that 
cannot be described in detail. Of course, the idea of a 
stochastic approach is hardly new 0, 0, 0] , but there 
are several different ways to implement it. The system we 
want to describe is the particle beam endowed with some 
measure of time-reversal invariance. For such a physical 
situation one is faced with a many-body dynamical sys- 
tem which is at the same time stochastic (due to the ran- 
dom collisions and the fast, irregular configurational dy- 
namics) and conservative (time-reversal invariant). One 
then needs a generalization of classical mechanics for con- 
servative deterministic systems to a classical mechanics 
of conservative stochastic systems, i.e. a stochastic me- 
chanics. 

Despite the rather widespread misconception that 
stochastic processes in dynamics are always associated 
to the description of dissipative and irreversible behav- 
iors, it is in fact possible for specific stochastic dynamical 
systems to exhibit time-reversal invariance. In general, 
these systems are such that a deterministic dynamics in 
terms of known external potentials is ascribed to a system 
whose kinematics, due to some intrinsic irregularities, is 
instead stochastic ^{|. This is the case we conside for 
the collisional beam. 

In previous papers [2(1 l2lj we analysed some basic 
properties of stochastic mechanics (SM) and we also in- 
troduced some basic criteria of mechanical stability in or- 
der to provide phenomenological support to the scheme 
of SM for classical stochastic dynamical systems. In par- 
ticular, in Ref. [2£j we applied SM to the description of 
the dynamics of charged particle beams, and we showed 
how the criteria of mechanical stability allow to connect 
the (transverse) emittance to the characteristic micro- 
scopic scales of the system and to the total number of 
particles in a given beam. Moreover, this method allows 
to implement techniques of active control for the dynam- 
ics of the beam. In Ref. [2(| these techniques have been 
proposed to improve the beam focusing and to indepen- 
dently change the frequency of the betatron oscillations. 

In the present paper, the SM of particle beams is used 
to investigate the nature, the size and the dynamical 
characteristics of beam halos. In particular, we deter- 
mine the effects of the space-charge interaction on the 
formation of the halo by calculating the growth of the 
beam emittance, and by estimating the probability to 
find particles far away from the beam core. We then 
consider the reverse problem i.e., by imposing reason- 
able forms of the stationary halo distributions we derive 
the effective spatial potentials that are associated to such 
distributions. 



The paper is organized as follows: in Sections |TTland lTTll 

we summarize the results of previous papers on the model 
of SM for particle beams with special emphasis on the 
techniques of active control needed to engineer a desired 
beam dynamics. We then apply SM to study the forma- 
tion of the beam halo. In Section lTVI we show how to treat 
with SM the effects of space charge due to the electro- 
magnetic interactions among the particles. We introduce 
and solve sclf-consistcntly a set of coupled, nonlinear hy- 
drodynamic equations for the beam density and phase, 
the space-charge potentials, and the external electromag- 
netic field. Solution of the coupled equations leads to the 
first main result of our paper, i.e. that the final trans- 
verse particle distribution is greatly broadened due to the 
presence of the space-charge potential. This broadening 
gives rise to a penetration of the distribution well into 
spatial regions very far from the core of the beam, thus 
realizing a type of continuous halo, i.e. without inter- 
mediate voids of particles (nodes of the transverse beam 
distribution). This effect is reflected in a corresponding 
growth of the transverse emittance that we estimate from 
the solution of the self-consistent equations. In Section 
[V]we study the reverse problem: given different station- 
ary transverse distributions of the halo around the core 
of the beam, such as a continuous broadened distribu- 
tion or a distribution with a gap (void) between the halo 
and the core, we determine the dynamics that leads to 
such configurations. This is important since it enables 
us to deduce the form of the effective potentials which 
control the stationarity of the halo distribution. Finally, 
in Section IVll wc draw our conclusions and discuss future 
directions of research. 



II. STOCHASTIC BEAM DYNAMICS 

The usual way to introduce a stochastic dynamical sys- 
tem is to modify the dynamics in phase space by adding 
a Wiener noise only in the equation for the momentum, 
in order to preserve the usual relations between position 
and velocity: 

mdQ(i) = P(t)dt 

dP(t) = F(Q(t),P(t),t)dt + /?dB(f) (1) 

where m is the mass of the system, Q(£) is the posi- 
tion variable, P(t) is the momentum variable, F is the 
external force, and B(t) is a fluctuating force modelled 
by a Wiener process; in general, this procedure yields a 
differentiable, but non Markovian process Q(i) for the 
position (here j3 denotes the autocorrelation of the fluc- 
tuating force). The standard example of this approach is 
that of a Brownian motion in a force field described by 
the Ornstein-Uhlenbeck system of stochastic differential 
equations (SDE) j^: 

mdQ(t) = P(t)dt 

dP(t) = K(Q(i), t) dt - m 7 P(t) dt + /3 dB(t) (2) 
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where K is the spatial gradient of the external poten- 
tial, and 7 is the damping constant. This system can be 
reduced to yield a non diffcrcntiable, Markovian process 
Q(i) acted on by a Wiener noise, with diffusion coeffi- 
cient D, to the equation for the position. In this case we 
are forced to drop the equation for the momentum since 
Q(t) is no longer diffcrcntiable, and we can reduce the 
stochastic system to the single SDE 

dQ(t) = v (+) (Q(t),t)dt + y/DdW(t) (3) 

where vr+> (r, t) is the so called forward velocity, and 
dW(t) = 'W(t+dt)—'W(t) is the increment of a standard 
Wiener noise W(i); as it is well known, this increment is 
a Gaussian process. The standard example of this reduc- 
tion is the Einstein- Smoluchowski approximation of the 
Ornstein-Uhlenbeck process in the overdamped case. 

It is however possible to introduce Eq. Q in a com- 
pletely different context, as the defining equation for a 
theory of stochastic dynamical systems in configuration 
space, rather than in phase space. This amounts to con- 
sider a system whose kinematics, rather than its dynam- 
ics, is taken ab initio to be random. By doing this, no 
external sources of dissipation and irreversibility are in- 
troduced on the forces acting on the system. Instead, a 
source of randomness is assumed that perturbs only the 
configurations of the system. For such systems with ran- 
dom kinematics, the dynamics can then be introduced 
either by generalizing the Newton equation |22l l23l, o r 
by introducing a stochastic variational principle |22l.l24j| . 
The crucial aspect is that such a dynamics turns out 
to be both stochastic and conservative. In this scheme, 
since Q(i) is not differentiable, there is no velocity as a 
standard derivative. Instead it is possible to define in a 
suitable way [2^ an average velocity v/_|_) in the forward 
time direction, and an average velocity in the backward 
time direction V(_) : these functions of r and t are differ- 
ent, and they both coincide with the usual velocity only 
if the process is differentiable, i.e. if the kinematics is 
deterministic. The relations between vv_|_) and vv_) will 
be introduced in the subsequent discussion. It is impor- 
tant to remark that vy_|_) (r, t) is no more an a priori given 
field: it plays now the role of a dynamical variable. 

The stochastic dynamical scheme sketched above is 
known as Stochastic Mechanics (SM), and most of its 
applications have concerned the problem of developing 
a classical stochastic model for the simulation of Quan- 
tum Mechanics. Nonetheless, it is a very general model 
which can be applied to very different stochastic dynam- 
ical systems endowed with time-reversal invariance |25j . 
We will show below how one can derive from the stochas- 
tic variational principle two coupled, nonlinear hydrody- 
namic equations for the density and the phase of a dy- 
namical system in configuration space. These two real, 
coupled equations can be recast into a single complex 
equation whose form is completely equivalent to that of 
a Schrodingcr differential equation. In this sense, some 
authors denote classical dynamical systems described by 
SM as quantum-like systems, in analogy with other re- 



cent studies on the collective dynamics of charged particl 
beams |2rll27j. 

In general, SM can be used to describe every conser- 
vative, stochastic dynamical system satisfying fairly gen- 
eral conditions: for instance, it is known that for any 
given conservative diffusion there is a correspondence be- 
tween the associated diffusion process and a solution of 
the Schrodingcr equation with Hamiltonians associated 
to suitable vector and scalar potentials . Under some 
regularity conditions this correspondence is one-to-one. 
The usual Schrodingcr equation, and hence quantum me- 
chanics, are recovered when the diffusion coefficient D 
is independent from the values taken by the conserva- 
tive diffusion process and coincides with h/2m, where 
h denotes the Planck constant. Here we are not inter- 
ested in the stochastic modelling of quantum mechanics, 
but rather in the classical stochastic dynamics of particle 
beams. In this instance the unit of action appearing in 
the effective Schrodinger equation is linked to the beam 
omittance as will be clarified below. 

We will now describe the stochastic process performed 
by a representative particle of the beam that oscillates, 
in a comoving reference frame, around the closed ideal 
orbit in a particle accelerator. We consider the three- 
dimensional (3d) diffusion process Q(t), taking the val- 
ues r, which describes the motion of the representative 
particle and whose probability density coincides with the 
particle density of the beam. The evolution of this pro- 
cess is ruled by the Ito SDE © where the diffusion co- 
efficient D is supposed to be constant. The quantity 
a = 2mD, which has the dimensions of an action, will be 
connected later to the characteristic transverse emittancc 
of the beam. Eq. @ defines the random kinematics per- 
formed by the collective degree of freedom, and replaces 
the simple deterministic kinematics 

dq(i) =v(q(f ),i)tft (4) 

for the differentiable trajectory q(f). 

We are in a situation in which we have both a ran- 
dom, diffusive kinematics and a time reversal invariance. 
Therefore, to counteract the dissipation, one must im- 
pose a conservative dynamics on the stochastic kinemat- 
ics, at variance with the purely dissipative Fokkcr-Planck 
or Langevin dynamics. A conservative dynamics imposed 
on a random kinematics can be introduced by a suitable 
stochastic generalization of the least action principle of 
classical mechanics [24[. This is achieved by replacing 
the classical deterministic kinematics (0J with the ran- 
dom diffusive kinematics © as the independent configu- 
rational variables of the Lagrangian density entering the 
action functional. The equations of motion thus obtained 
by minimizing the action functional take the form of two 
coupled hydrodynamic equations that describe the dy- 
namical evolution of the beam profile and of the associ- 
ated velocity field. In the following we briefly sketch the 
derivation of the hydrodynamic equations, referring for 
details to references [22, |24| . 

Given the SDE © , we consider the probability density 
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function (pdf) p(r,t) associated to the diffusion Q(i) so 
that, besides the forward velocity V(_|_)(r, t), we can also 
define the backward velocity 

V( _ ) (r ) t)^v (+) (r,i)-2i?^M. (5 ) 

It is also useful to introduce the current and the osmotic 
velocity fields, defined as: 



_ v (+) 



'(-) 



v (+ ) - v r _ 



(-) 



The velocities in Eq. © have a transparent physical 
meaning: the current velocity v represents the global 
velocity of the density profile, being the stochastic gen- 
eralization of the velocity field of a classical perfect fluid. 
On the other hand the osmotic velocity u is clearly of 
intrinsic stochastic nature, for it is a measure of the non 
differentiability of the stochastic trajectories, and it is 
related to the spatial variations of the density. In the 
limit of a deterministic process, i.e. of a diffusion that 
tends to a deterministic, differentiable trajectory q(i), 
the current velocity v(r, t) tends to the classical veloc- 
ity field v(q(i),t), and the osmotic velocity u tends to 
zero (the trajectory becomes differentiable, so that the 
forward and backward velocities, i.e. the left and right 
derivatives coincide). 

In order to establish the stochastic generalization of 
the least action principle one introduces a mean classi- 
cal action averaged over the probability density function 
of the diffusion process, in strict analogy to the classi- 
cal deterministic action. In fact, the main difficulty in 
the stochastic case is due to the non differentiable char- 
acter of the sample paths of a diffusion process which 
does not allow to define the time derivative of the pro- 
cess Q(i). Hence the definition of a Lagrangian density 
and of an action functional is possible only in an aver- 
age sense through a suitable limit on expectations. The 
stochastic action is then defined as 123, |2J1 



A 



to 



lim E 



dt, (7) 



where E( • ) = J( ■ )p(r, t) dr denotes the expectation of a 
function of the diffusion process Q(t), V is an external 
potential, and AQ(f) = Q(t + At) - Q(t) is the finite 
increment of the process. It can be shown that the mean 
action (JJJ associated to the diffusive kinematics (0 can 
be recast in the following particularly appealing Eulcrian 
hydrodynamic form j2^| : 



A 



dt / dr 



(v 2 



V 



P(r,*), (8) 



where v and u are defined in Eq. © • This Eulcrian form 
of the action immediately shows that when u = (and 
hence the trajectories are differentiable) it coincides with 



the action for the conservative dynamics of a classical Li- 
ouvillc fluid. The stochastic variational principle now fol- 
lows by imposing the stationarity of the stochastic action 
(6 A = 0) under smooth and independent variations dp of 
the density, and (5v of the current velocity, with vanishing 
boundary conditions at the initial and final times. 

As a first consequence we get that the current velocity 
has the following gradient form: 



mv(r,t) = VS(r,t) , 



(9) 



which can be also taken as the definition of the phase S. 
The non linearly coupled Lagrange equations of motion 
for the density p and for the current velocity v, of the 
form ©, are a continuity equation typically associated 
to every diffusion process 

d t p - -V • (pv) , (10) 

and an evolution equation for the conservative dynamics 



d t S 



2mD 



■ y 2 yp 

VP 



+ V(r,t) =0. 



(11) 



The above equations give a complete characterization 
of time-reversal invariant diffusion processes. The last 
equation is formally analogous to the Hamilton-Jacobi- 
Madclung (HJM) equation originally introduced in the 
hydrodynamic formulation of quantum mechanics by 
Madelung [2^. However, the physical origins of the two 
equations are profoundly different, as in the quantum 
mechanical case the diffusion coefficient is related to the 
fundamental Planck constant by the relation D = h/2m. 
Notice that when D = (namely when the kinematics is 
not diffusive) the equation i|ll|) coincides with the usual 
Hamilton-Jacobi equation for a Liouvillc fluid. Since © 
holds, the two equations i|10|) and (|11|) can be recast in 
the following form 



dtp 



1 



V ■ ( P VS) 



m 
2m 



(12) 



2 m D 2 ^-^ - V(r,t) (13) 



which now constitutes a coupled, non linear system of 
partial differential equations for the couple (p, S) which 
completely determines the state of our beam. On the 
other hand, because of JSJ, this state is equivalently given 
by the couple (p, v). 

Eqns. i|12|) and i|13|) describe the collective behavior of 
the beam at each instant of time through the evolution 
of both its particle density and its velocity field. It is 
important to notice that, introducing the representation 



with 



a = 2mD , 



(14) 



(15) 
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the two coupled real equations (|12f> and (|13fl are equiv- 
alent to a single complex Schrodinger equation for the 
function ip, with the Planck action constant H replaced 
by the unit of action a: 



iadtip 



2m 



V 2 tp + Vip. 



(16) 



In this formulation the phenomenological "wave func- 
tion" ip carries the information on the dynamics of both 
the beam density and the beam velocity field, since the 
velocity field is determined via Eq. @ by the phase func- 
tion S(r,t). This shows, as previously claimed, that our 
procedure, starting from a different point of view, leads 
to a description formally analogous to that of the so- 
called quantum-like approaches to beam dynamics |26j| . 
We remark again that obviously Eq. (|16fl has not the 
same meaning as in quantum mechanics. The parameter 
a defined in Eq. 1)15(1 has the dimensions of an action, but 
in general a ^= h. In fact, a is not a universal constant 
and its value cannot coincide with h, as its value depends 
on the physical system under study. However, a plays in 
some sense a role similar to that of h in quantum me- 
chanics since, as we will see in the following, it identifies 
a lower bound for the beam omittance in phase space, 
and gives a measure of the position-momentum uncer- 
tainty arising from the stochastic dynamics of the beam. 
Thus the Schrodinger equation l|16(l for the stochastic 
mechanics of particle beams presents some features rem- 
iniscent of quantum mechanics, but at the same time is 
a deeply different theory with a different physical mean- 
ing. In particular, while in quantum mechanics a system 
of TV particles is described by a wave function in a 3N- 
dimensional configuration space, in the scheme of SM 
the normalized probability density distribution ^(r)) 2 is 
a function of only the three space coordinates in phys- 
ical space. If there are N particles in the beam, then 
the function N\ip(r)\ 2 is simply the particle density in 
physical space. 



III. CONTROLLED BEAM STATES 

In the previous section we have introduced two cou- 
pled equations that describe the dynamical behavior of 
the beam: the first is the ltd equation J3J, or cquiva- 
lently the continuity equation (cither (|10|l or (|12|l 'l: the 
second is the Hamilton-Jacobi-Madelung (HJM) equa- 
tion (either (|11|) or flU} ). Here we briefly summarize a 
general procedure, already exploited in Refs. [20ll3Gll3l| . 
to engineer a controlled dynamics of stochastic systems. 
We will then apply this method to the control of the 
beam halo in Section Ivl 

To this end it is useful to show, by simple substitution 
from © , that Eq. (|10fl is equivalent to the Fokker-Planck 
(FP) equation 



d t p= - V- [v (+) p] +DV 2 p, 



(17) 



formally associated to the Ito equation ©. The HJM 
equation (|ll|l can also be cast in a form based on V( + ) 
rather than on v, namely: 



d t S 



mDv 



-mD 2 



(+)■ 



V 2 p 



Vp 



Vp 
P 



V 



y v (+) 



mD V( + ) V In / 
+mD 2 V 2 In / - V , 



(18) 



where we have defined the function / = p/JV , and Af is a 
constant such that / is dimensionless. On the other hand, 
we know from Eqns. (0) and © that also the forward 
velocity V( + ) is irrotational: 



v (+ )(r,i)=V%i), 



(19) 



where the scalar functions W and S are connected by the 
relation 

S(r, t) = mW(r, t) - mD In /(r, t) - 6(t), (20) 

and 9 is an arbitrary time-dependent function. Through 
Eq. I|19|) . Eqns. I|17|) and l|18|) can be cast in a form anal- 
ogous to the system of Eqns. H12|) and 1|13|) . In this case, 
the couple of functions which determine the state of the 
beam is (p, W) (or cquivalently (p, V( + ))). 

It is worth noticing that the possibility of a time rever- 
sal invariance [22j is assured by the fact that the forward 
velocity V( + ) (r, t) is not an a priori given field, as is the 
case for dissipative diffusion processes of the Langevin 
type. Rather, it is dynamically determined at any in- 
stant of time, for any assigned initial condition, through 
the HJM evolution equation (|11|) . 

Let us suppose now that the beam density p(r, t) is 
given at some time t. We may think for instance of an 
engineered evolution from some initial density toward a 
final, required state with suitable properties. In particu- 
lar, we could imagine to have a halo- ridden particle dis- 
tribution in the beam that should be steered toward a 
final halo-free distribution. We know that p must be a 
solution of the FP equation Ijl7[l for some forward veloc- 
ity field V( + )(r, t), which we consider here as not given a 
priori. Since also Eq. (|19fl must be satisfied, this means 
that - for a given p - we should find an irrotational V( + ) 
in such a way that the FP equation (|17fl is satisfied. In 
other words, we are required to solve the partial differ- 
ential equation 



V-v 



(+) 



Vp 



'(+) 



= D 



V 2 p d t p 



— , (21) 
P 



for the irrotational vector field V( + ) , or cquivalently, by 
taking Eq. (|19fl into account, the second order partial 
differential equation for the scalar field W: 



\7 2 W + — P - ■ VW 



D 



V 2 p dtp 



(22) 
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This is not, in general, an easy task, but we will give here 
the solutions for a few simple cases that will be useful in 
the following. 

In the one-dimensional (Id) case (p{x, t)), the FP equa- 
tion defined for a < x < b (here we could possibly have 
a = — oo, and b = +00) is: 



with 



d t p = -d x [v {+) p] + Dd 2 x p = d x J, 



(23) 



where J = Dd x p — vt + \p is the probability current. It is 
easy to see now that the solution is 



v {+) (x,t) = - 



c(t)+Dd x p- / d t p(x',t)dx' 



(24) 



where c(t) is an arbitrary function. Moreover, since the 
conservation of probability imposes, in particular, that 
J(a) = 0, it is easy to sec that we must choose c(t) = 
so that finally we have 



v (+)( x >t) = D 



d x p _ 1 
P P 

Of course, in this case we also have 



dtp(x',t) dx' . 



W(x,t) = / v {+) {x' ,t)dx' , 



(25) 



(26) 



This Id solution will be useful later for a simplified study 
of beam distributions when we take x as one of the trans- 
verse coordinates of the beam. 

Another case of practical interest is that of a 3d sys- 
tem with cylindrical symmetry around the z-axis (which 
will be supposed to be the beam longitudinal axis). If we 
denote with (r, ip, z) the cylindrical coordinates, where 
r = \J x 2 + y 2 , we will suppose that p(r, t) depends only 
on r and t, and that V( + ) = V(+-\ (r, t) r is radially di- 
rected, with modulus depending only on r and t. In this 
case, it is straightforward to see that the equation for 
v (+) is 



d. 



9 r P 



whose solution is 

«(+)( r > *) 
In this case 



D 



d r p 



1 



D 



d 2 P | d r p 

P rp 



d t p 
P 



r P Jo 



d t p(r',t)r' dr' . 



(27) 
(28) 

(29) 



Finally, in the 3d stationary case the beam density p is 
independent of t. This greatly simplifies the treatment, 
as it is easy to check by simple substitution, and we sim- 
ply find 



W(r,t) = / v {+) (r',t)dr' , 



v ( +)(r)=D 



Vp(r) 
p(r) ' 



(30) 



W(r) =Dln/(r) 



(31) 



This result corresponds to the well known fact that for 
stationary forward velocities the FP equation (|17|l always 
admits stationary solutions of the form 



(32) 



Remark however that also for particular choices of a 
non stationary p(r, t) it is sometimes possible to find a 
stationary velocity field V( + )(r) such that the FP equa- 
tion (|17|l is satisfied (think for instance to the Ornstein- 
Uhlcnbeck non stationary solutions): this fact will be of 
use in the following Sections. 

Now that p and W (or cquivalently V( + )) are given and 
satisfy Eq. JT7J|, we should also remember that they will 
qualify as a good description of our system only if they 
also are a solution of the dynamical problem, namely if 
they satisfy the HJM equation (|18|) . Since W, vc+) and 
p are now fixed, this last equation must be considered 
as a constraint relation defining an external controlling 
potential V which, after straightforward calculations and 
in terms of the dimcnsionlcss distribution /, turns out to 
be of the general form: 

V(r,t) = ml) 2 V 2 In/ + mL>(<9 t hi/ + v (+) . Vln/) 



--v 2 +) -md t w + e. 



(33) 



This procedure can in principle be applied also to more 
complicated instances, for example to engineer a beam 
dynamics that keeps the beam coherent even in the pres- 
ence of aberrations, or to engineer some desired evolu- 
tion. However the solutions will not always be available 
in closed form, and will in general require some approx- 
imate treatment. Here we present the simple analytic 
solutions for the three particular cases discussed above. 

In the Id case f(x, t), with given by Eq. (|25|l . the 
controlling potential is: 

V(x,t) = 6 + m.D 2 dl In / + mD (d t In f + v {+) d x In/) 

m i x 
-jvf +) -mJ d t v i+) (x',t)dx'. (34) 

In the 3d, cylindrically symmetric case f(r,t), with vr+) 
given by (|28|l . we have instead: 



V(r,t) 



d r (rd r In /) + mD(d t In / + «( + )<9 r In /) 



-rrv 2 {+) -m / d t v (+) (r',t)dr' + 6. (35) 



Finally, in the 3d stationary case p(r) with V( + ) given by 
Eq. (|30|l . the potential is still time-dependent because of 
the presence of the arbitrary function 9(t), and it reads: 



V(r,t) = 6 + 2mD' 



\fP 



(36) 
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In order to make it stationary it will be enough to require 
that 6 it) be a constant E, so that 



VYr = E + 2mD 2 — ±Jl 



(37) 



In this stationary case the constant E fixes the zero of the 
potential energy, and the phcnomcnological wave func- 
tion (|14|l takes the form 

V(r,t) = yfpe- iEt ' a (38) 
typical of stationary states. 

IV. SELF-CONSISTENT EQUATIONS 

One of the possible mechanisms for the formation of 
the halo in particle beams is that due to the unavoid- 
able presence of space charge effects. In this Section 
we will investigate this possibility in the framework of 
our hydrodynamic-stochastic model of beam dynamics. 
To this end, we take into account the space charge ef- 
fects by coupling the hydrodynamic equations of stochas- 
tic mechanics with the Maxwell equations which de- 
scribe the mutual electromagnetic interactions between 
the particles of the beam. We thus obtain a self- 
consistent, stochastic magnctohydrodynamic system of 
coupled nonlinear differential equations that can be nu- 
merically solved to show the effect of the space charge. 

In the following, the reference physical system will be 
an ensemble of TV identical copies of a single charged par- 
ticle embedded in a particle beam and subject to both 
an external and a space-charge potential. In a reference 
frame comoving with the beam, our system is then de- 
scribed by the Schrodinger equation l|43|l . where a is the 
unit of action (emittance) and H the Hamiltonian oper- 
ator which will be explicitly determined in the following. 
Since in general ip is not normalized, we introduce the 
following notation for its constant norm 



Ml 



R 1 



|V(r,t)| 2 d J r, 



(39) 



so that, if N is the number of particles with individual 
charge qo, the space charge density of the beam will be 



p ac (r,t) = Nq 



(40) 



On the other hand its electrical current density will be 

faM.^ isffl^M , (4i) 

m 



which vanishes when the wave function is stationary, 
namely when 



ip(r,t) =uir)e- lEt ' a 



(42) 



For a charged particle in the beam the electromagnetic 
field is the superposition of the space-charge potential 
(A sc , $ sc ) due to the presence of p sc and j sc , plus the ex- 
ternal potential (A ea;t , <& ex t), and hence the Schrodinger 
equation takes the form 



iadt/ip = 



1 



2mc 2 



[iacV - go(A sc + A ext )} ip 
-h7o(3 > sc + 3w)V' ■ 



(43) 



Eq. I)43|) has then to be coupled with the Maxwell equa- 
tions for the vector and scalar potentials (A SC ,<E> SC ). 
Since we are in a reference frame comoving with the 
beam, we can always assume that the wave function is of 
the stationary form (|42|l . In this case we have j sc = 0, 
so that A sc = 0. Finally, by taking also A ext = 0, our 
system is reduced to two coupled, nonlinear equations 
for the pair (it, $ sc ). For details see Appendix^] If 
the beam with space-charge interactions stays cylindri- 
cally symmetric the function u will depend only on the 
cylindrical radius r and our system becomes (see equa- 
tions HA5|1 and HA6f) in Appendix IaT) : 



Eu 



a 

2m 



+ — + iy ext + v sc )u, 



2ne LA 



u 2 = 



= - V. 



V 

v sc 

r 



where 



A = I ru (r) dr , 
lo 



2itLA. 



(44) 
(45) 

(46) 



with L the length of the beam, and 

V ext (r,t) = q $ ext (r,t) , V sc (r, t) = q $ sc ir, t) . (47) 

Eq. (|44fl is the stationary Schrodinger equation in the ex- 
ternal scalar potential, and Eq. I|45|) is the Poisson equa- 
tion for the space-charge scalar potential. We proceed 
to solve numerically this nonlinear system and compare 
its solution with that for a purely external potential V ex t 
which is a cylindrically symmetric, harmonic potential 
with a proper frequency u in absence of space charge (see 
Appendix |B|) : V ext = (mw 2 r 2 )/2, with r = \J x 2 + y 2 . 
Let us first introduce the dimensionlcss quantities 



*V2' 



wis) 



vis) 



w 



aV2 



aV2 



= cr 3/2 w(r) 



(48) 



Ama' 



Vsdr) , 



where a 2 = a/2mu> is the variance of the ground state of 
the cylindrical harmonic oscillator without space charge. 
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Eqns. H44f) and l|45|l can then be recast in the dimension- 
less form 

sw"(s) + w'{s) + \/3 ~ s 2 - v(s)} sw{s) = , (49) 



sv"(s) + v'(s) + b sw 2 (s) = , 



(50) 



where 



p = 


4ma 2 
o 2 


2E 
E=—, 




b = 


4m<7 2 


Nq 2 1 _ 2 


Nq 2 1 


n 2 


27reoi B au 


27re £ B 


B = 


Aa 
2 


poo 

/ sw 2 (s)ds. 
Jo 





o=— 5Z7To. (51) 



The effect of the space charge will be finally accounted 
for by comparing the normalized solution w(s) with the 
unperturbed ground state of the cylindrical harmonic os- 
cillator (see Appendix lB|l : 



,-r 2 /4cr 2 



(52) 



We have solved numerically the system (|49ll , (|5UII by ten- 
tatively fixing one of the two free parameters b and ft, and 
then searching by an iterative trial and error method a 
value of the other such that the solution shows, in a given 
interval of values of s, the correct infinitesimal asymp- 
totic behavior for large values of s. We have then nor- 
malized the solutions w(s) by calculating numerically the 
value of B. It is clear from the definition of the dimcn- 
sionless parameters B, b and j3 l|51|) that the value of f} is 
a sort of reduced energy eigenvalue of the system, while 
the product 



l = Bb, 



(53) 



which is by definition a non negative number, will play 
the role of the interaction strength, since it depends on 
the space-charge density along the linear extension of 
the beam. Reverting to dimensional quantities, since 
a 2 /2ma 2 has the dimensions of an energy, the two rele- 
vant parameters are 



E = (3 



4toct 2 



Nql 
2t:€qL 
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Ama 2 



(54) 



which are respectively the energy of the individual parti- 
cle embedded in the beam, and the strength of the space- 
charge interaction. 

In the following we will limit ourselves to discuss so- 
lutions of Eqns. 149|) . and (|5L)fl without nodes (a sort of 
ground state for the system). It is possible to see that 
no solution without nodes can be found for values of the 
space-charge strength 7 beyond about 22.5 and that for 
values of 7 ranging from to 22.5, the energy j3 decreases 



2 

1.5 
1 

0.5 



10 15 
9 



20 



FIG. 1: The energy /3 as a function of the space-charge 
strength 7. Actually these quantities are dimensionless, mea- 
sured in units of the energy a 2 /4ma 2 ~ 37.5 eV. 



monotonically from 2.0000 to -0.0894. If the unit of ac- 
tion a is fixed at a given value, it is apparent that the 
value of 7 is directly proportional to the charge per unit 
length Nqo/L of the beam: a small value of 7 means a 
rarefied beam; a large value of 7 indicates that the beam 
is intense. A halo is supposed to be present in intense 
beams, while in rarefied beams the behavior of every sin- 
gle particle tends to be affected only by the external har- 
monic potential. 

It is useful to provide a numerical estimate (in MKS 
units) of the relevant parameters. We are assuming a 
beam made of protons, so that m and qo are the pro- 
ton mass and charge, while eo is the vacuum permittiv- 
ity. The parameter a is determined by the strength of 
the external potential, and it is a measure of the trans- 
verse size of the beam when no space charge is taken 
into account, so that the ground state has the unper- 
turbed form Ij52(l . From empirical data, a reasonable es- 
timate yields a w 10~ 3 m. On the other hand the value 
of N/L clearly depends on the particular beam we are 
considering: usually a value of 10 11 particles per meter 
is considered realistic. As for the parameter a (|f 5J) we 
have already suggested in Section [I] that its value can 
be connected to the beam emittance. On the ground 
of accepted experimental values of the beam emittance 
(usually measured in units of length) we can assume 
a/mc of the order of 10 m. In fact, taking m, qo, a, 
and N/L fixed at the above-mentioned values, and as- 
suming a space-charge strength 7 moderately large, i.e. 
2.0 < 7 < 22.5, we can determine a from (|54|) and we 
get 7.2 x 10~ 7 < a/mc < 1.3 x 10~ 7 , a number which 
is in good agreement with the experimentally measured 
values of the beam emittance. From now on we will fix 
a at its approximate central value: 



— w 4.0 x 10" 

mc 



(55) 



Since m and qo have the proton values, and a « 10~ 3 m, 
then the assigned values of the dimensionless parameters 
7 and (3 will fix respectively, the values of N/L and E by 
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s w 2 HsL 




s w 2 HsL 



FIG. 2: The radial distribution sw 2 (s) compared with the 
distribution in absence of space charge (dashed line). The 
space-charge strength is 7 ~ 11.2. All quantities are dimen- 
sionless. 



FIG. 4: The radial distribution sw 2 (s) compared with the 
distribution in absence of space charge (dashed line). The 
space-charge strength is 7 ~ 22.5. All quantities are dimen- 
sionless. 





FIG. 3: 3-di mensional view of the transverse distribution 
w 2 (^/x 2 + y 2 ) for a space-charge strength of 7 ~ 11.2. 



FIG. 5: 3-di mensional view of the transverse distribution 
w 2 {\J x 2 + y 2 ) for a space-charge strength of 7 ~ 22.5. 



Eqns. (|54() . This implies that by changing the beam in- 
tensity (namely N/L and 7) one correspondingly changes 
the energy (namely E and (3) of the individual particle 
embedded in the beam. In Fig. ^ we show the behavior 
of f3 as a function of 7. The quantities (3 and 7 are di- 
mcnsionlcss: the true physical quantities (energies) are 
obtained by multiplicating them by the unit of energy 



a 



37.5 eV. 



(56) 



The overall effect of the space charge in this model is 
a conspicuous spreading of the transverse distribution of 
the particles in the beam with respect to the unperturbed 
ground state distribution When 7 = the poten- 

tial due to space charge vanishes, and the solution ex- 
actly coincides with the ground state of the cylindrical 
harmonic oscillator with variance a 2 . When 7 > the 
transverse distribution begins to spread, as we show in 



Figrs. and for a value of the space-charge strength 
7 w 11.2, and in Figrs. 0] and [5] for a value of the space- 
charge strength 7 ss 22.5. Notice that when comparing 
distributions, the true radial (dimcnsionlcss) density is 
s w 2 (s) and not just w 2 (s). Recall that w(s) is a normal- 
ized solution of Eqns. (gHJ) and (JSUJl. In Figrs. El and [7| 
we show the radial transverse form of the total poten- 
tial, i.e. external plus space-charge given by the solution 
v(s) of Eqs. (|4"9|) and (j5TJ|) . for the values 7 = 11.2 and 
7 = 22.5 of the space-charge strength. To give a more 
quantitative measure of the flattening and broadening of 
the transverse distribution we compare numerically the 
probabilities of finding a particle at a relatively large dis- 
tance from the beam longitudinal axis with and without 
space charge. The quantity 



P 7 (c) = 



/V2 



s w 2 (s) ds 



(57) 
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FIG. 6: The total radial potential felt by a particle in the 
beam compared with the harmonic potential (dashed line) 
for a space charge strength of 7 w 11.2. 




FIG. 7: The total radial potential felt by a particle in the 
beam compared with the harmonic potential (dashed line) 
for a space charge strength of 7 « 22.5. 



is the probability of finding a particle at a distance 
greater than ca for systems with a given strength 7 of the 
space charge coupling. For instance, considering the two 
different situations 7 = (no space charge) and 7 = 22.5 
(strong space charge) we have: 



P (10) ps 1.9 x 10" 



22.5 



(10) « 1.7 x 10" 6 . (58) 



We see that the probability of finding particles at a dis- 
tance larger than 10cr from the core of the beam is en- 
hanced by space charge by many orders of magnitude. 
This means, for example, that if in the beam there are 
10 11 particles per meter, while practically no one is found 
beyond 10er in absence of space charge, for very strong 
space-charge intensity we can find up to 10 5 particles per 
meter at that distance from the core. 

The above analysis shows that in the hydrodynamic- 
stochastic theory of charged beams, the space-charge po- 
tential induces a strong broadening of the unperturbed 
transverse density distribution of the beam, thus yield- 
ing a small, but finite probability of having particles at 
a distance well away from the core of the beam. 



The hydrodynamic-stochastic model allows also an es- 
timate on the growth of the omittance due to the presence 
of the space charge. The emittance can be calculated by 
exploiting a structure of uncertainty products that is in- 
herent to the SM. In particular, the transverse emittance 
can be calculated as Ax • Ap x , where we have: 



Ax = V(z 2 ) - (x) 2 , A Px = V(p2) - (p x )< 



(59) 



More explicitly: if w(s) is a solution of the coupled dy- 
namical equations normalized in the sense that 



s w 2 (s) ds = 1 



the corresponding 3d probability density is 

p(r,<p,z) = u 2 (r)H - \z\ 

where H(s) is the Heaviside function and 

1 



(60) 



(61) 



u{r) 



This probability density function is normalized since 



(62) 



+ OO p27T 

r dr I dip 



dz p(r, ip, z) = 1 , 
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It is easy to show that the first two moments are: 



0. 



- / s 3 w 2 (s)ds, 
'0 



(63) 



(64) 



so that the root mean square deviation for the position 
is: 



Ax = \J {x 2 ) = a ]j~J s 3 w 2 (s) ds 



(65) 



As for the momentum p, we recall that in SM one can- 
not introduce a probability distribution directly in phase 
space. The momentum can nevertheless be recovered 
from the velocity field which is well defined in SM. Since 
we are considering a stationary state, only the osmotic 
part of the velocity field is non zero. Then the momen- 
tum field reads 



p(r) = a 



Vp(r) 



2a 



Vu(r) 



(66) 



p(r) u(r) 
As a consequence, the x-component of the momentum is 



Px = 2a 



d x u 



2a 



x u'(r) 



2a cos tp 



u'{r) 



u r u(r) u{r) 

and the first and second moments read 



(p*> =0, 



(pi) = — w' 2 (s)sds. 
a Jo 



(67) 



(68) 
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FIG. 8: The transverse dimensionless emittance e vs. the 
space-charge strength 7. The true dimensional emittance, 
calculated as position-momentum uncertainty product, is just 
£ = eat. 



Hence the root mean square deviation for the momentum 



is: 




„/2 



(s)s ds 



(69) 



We can then define the emittance £ as the position- 
momentum uncertainty product in phase space in the 
following way: 



£ 




s 3 w 2 (s)ds- / w' 2 (s)sds. 



(70) 

The emittance (|70[) can be thus estimated from knowl- 
edge of the numerical solutions w(s) and their first 
derivatives w'(s). It can be seen that the emittance is 
exactly a for 7 = (namely in absence of space charge), 
and grows with 7 to a value « 1.2 x a for 7 ss 22.5. This 
result is consistent with the expected growth of emittance 
produced by space-charge effects. The dimensionless ra- 
tio e = £ I a = Ax ■ Ap x /a is plotted as a function of 
7 in Fig. [5] This figure and the above discussion pro- 
vide evidence that in the model of SM of charged beams 
the constant a plays the role of a lower bound for the 
phase-space emittance. 



we should try to analyze the solutions with nodes of our 
radial self-consistent equations. Being this a time con- 
suming task, we resort for the time being to a different, 
simplified approach to study the structure and properties 
of beams of the ringed type. 

Moreover, the scheme introduced in the previous Sec- 
tion, although in principle very general, has the drawback 
that it does not allow for analytical expressions of veloc- 
ities and potentials, and therefore cannot suggest the en- 
gineering of controlling external potentials to remove the 
halo. To this end, in this Section we will develop a dif- 
ferent approach, still based on SM, by considering some 
preassigned density distributions for a beam with halo, 
and then looking for the kind of dynamics that produces 
them. 

Firs of all we will determine the analytical forms w(s) 
of the dimensionless radial distribution that can best ap- 
proximate the numerical, self-consistent solutions w(s) 
found in the previous Section. We will do that by choos- 
ing a family of trial functions w(s) and by subsequently 
minimizing the mean square error (m.s.e.) with respect 
to a given solution w(s). Since the radial component of 
the cylindrically symmetric (I = 0) solutions of the har- 
monic oscillator (|Blf) contain only even powers of the ra- 
dial coordinate (sec Appendix B), wc take as normalized 
trial functions 



w(s) 



,-s 2 /2cr 2 



{l + as 2 + bs 4 ) 



(71) 



with cr, a and b as free parameters, and 
e- s2 /° 2 (l + as 2 + bs 4 ) 2 sds 



aa 4 + (a 2 + 26)cr 6 + 6aba s + 12&V (72) 



as normalization constant. We then take one numerical 
solution w(s) of l|49|l and l|50|) for a given value of 7 and 
consider the mean square error 



J 



\w(s) — w(s)\ 2 s ds . 



(73) 



and numerically minimize its value by varying the pa- 
rameters cr, a and b. As an example, for 7 ss 19.3, we get 



V. STATIONARY HALO DISTRIBUTIONS 

The coupled nonlinear equations introduced in the pre- 
vious Section allow to introduce and explain a possi- 
ble mechanism of halo formation due to space-charge ef- 
fects. This mechanism, when considering solutions with- 
out nodes, amounts to the broadening of the transverse 
density distribution of the beam away from the beam 
core. This is a type of "continuous" halo. If we want 
to investigate the formation and the properties of other 
possible forms of halos - for example ring-like halos sep- 
arated from the core by an almost void space region - 



a « 1.81; aw 0.28; b « 0.01 ; c w 5.09 . (74) 

The corresponding density distribution (solid line) is 
compared in Fig.[!|]with the density distribution obtained 
by numerically solving the system Ij49(l and l|50(l (dashed 
line). From this explicit form of the approximating w(s) 
we can now also calculate the expression of the station- 
ary potential which controls the state. From Eq. I|37|l in 
the case of a cylindrically symmetric approximating state 
with 

u(r) = ^W), (75) 
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FIG. 9: The best approximation in mean square (solid line) of 
a numerical solution for the self-consistent equations (dashed 
line) for trial functions of the form 171H with 7 « 19.3. 



we have the controlling potential 



a 

2m 



u"{r) t lu'jr) 
r u{r) 



u(r) 



V(r) =E 

Introducing the usual dimensionless quantities 

w(s) = a 3/2 u(r) , 



(76) 



~ Ama 2 ~ 



(77) 



we may define the dimensionless controlling potential 

4mcr 2 - 



v(s) 



+ 



0-V2/ a 

w"(s) + l w'(s) 
w(s) s w(s) 



V(r) 



(78) 



For an approximating amplitude of the form Q71[l. and 
with the same values (I74() of the parameters, the control- 
ling potential reads 

, 0.442 - 0.629s 2 + 0.062s 4 - 0.001s 6 

v(s) = = . (79) 

v ' 1 + 0.284s 2 - 0.011s 4 v ' 

This potential is shown (solid line) in Fig. ^| where it 



is compared with the potential 



v[s) (dashed line) 



solution of Eqns. J3j5J| and JSUJ). 

The analytical approximations show a small, but rel- 
evant, difference with respect to the numerical solutions 
of Eqns. (|4*9")l and (J5DJ: w(s) has a node at about s « 5.4, 
and correspondingly the potential v(s) has a singularity 
at that point (out of the s range of Fig. llO|) . Instead, the 
numerical solutions w(s) and v(s) of Eqns. I|49|) and l|50|) 
show no such behavior. For values of s beyond s w 5.4, 
w(s) shows a small bump, not visible in Fig. because 
of the large scale. We show a zoom-up of it in Fig. ^2 
where the node is now clearly visible. This is a feature 
that we will find as well in other different approaches to 
the description of the beam distribution, and that could 
be tentatively connected with the presence of a halo spa- 
tially separated from the core of the beam. 



FIG. 10: The control potential (solid line) obtained from the 
best mean square approximation w(s) to a numerical solution 
for the self-consistent equations 1491 and 150H . for trial func- 
tions of the form 1711 1 with 7 = 19.3. It is compared with 
the potential s + v(s) (dashed line) obtained as numerical 
solution of and (5J5J). 



s w HsL 




FIG. 11: Blown up picture of Fig. EJthat shows the existence 
of a node in the position density distribution at s ~ 5.4. 



In our quest for an analytic description of the halo 
distribution we now depart from the scheme of the self- 
consistent nonlinear equations I49|) and JSUJ. Instead, 
while still assuming the space charge as the dominant 
factor, we introduce a simplified analytical model for the 
halo production. As a first simplification, valid for small 
values of the space-charge strength 7, we assume a frozen 
core for the beam. A frozen core is such that it produces 
a space-charge effect, induces a spreading of the density 
distribution, and is left unchanged by this spreading. In 
this case the potential produced by the charge distribu- 
tion is a priori given by the distribution of the frozen 
core and it simply enters linearly in the phenomcnolog- 
ical Schrodinger equation 1|16|) . As a consequence, the 
latter is now decoupled from the Maxwell equations for 
the electromagnetic field. The frozen charge distribution 
will be assumed cylindrically symmetric and transver- 
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sally Gaussian: 



Q(r) = Nq - 



-(x 2 +y 2 )/2a 



27TCT 2 



(80) 



where N is the number of particles with elementary 
charge qo, L is the longitudinal extension of the beam, 
and H(-) denotes the Heavisidc function. The cylindri- 
cal symmetry allows to determine via Gauss theorem the 
potential energy V sc (r) of a test particle of charge qo em- 
bedded in this charge distribution. It is a simple exercise 
to show that 



V sc (r) 
where C R 



2 r 



Ei - 



In- 



2ire L [ ' \ 2a 2 J 2a 2 
0.577 is the Euler constant, and 



Ei(x) 



t 



dt 



(81) 



(82) 



is the exponential-integral function. If we suppose that 
the test particle is acted upon by this space-charge poten- 
tial, and by the external, cylindrical harmonic oscillator 
potential (|B1() - which, without space charge, would keep 
the particle in the Gaussian state with variance a 2 - the 
total potential to be inserted in the Schrodingcr equation 
(Unj reads 



V(r) 



8ma 4 



Nq 2 p 

2T76QL 



Ei - 



2a 5 



-In 



2cr 2 



(83) 



In terms of the usual dimensionless quantities we can also 
write 



Ama 2 Nq { 



ay/2' 



t 2 2tt£[)L 



v(s) 



Ama 



,oV2/ 
s 2 +7[Ei(-s 2 ) 



V(r) 



Ins 2 



A plot of the dimensionless potential v(s) for 7 
shown in Fig. ^| 

For stationary solutions of the form 



ip(r,t) =u(v)e- lEt/c 



Eq. 1|16|) becomes 



Eu{v) = Hu(r) = - — V 2 w(r) + V(r) u(r) , 
2m 



(84) 
3 is 

(85) 
(86) 



and due to the cylindrical symmetry of V it reduces to 
■ 2 r u'{r)~ 



Eu{r) = 

2m 



u"{r) 



+ V{r)u(r). (87) 
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FIG. 12: The dimensionless space-charge potential v(s) pro- 
duced by a frozen, cylindrically symmetric beam distribution 
with 7 = 3. 



Considering the dimensionless quantities 
Ama 2 



E , w(s) = w 



a' \aV2 
we finally obtain the reduced equation 



a 3/2 u(r) (88) 



(3w(s) = -w"(s) 



w'{s) 



+ v{s)w{s). (89) 



Despite its simple graphical behavior displayed in Fig.lT21 
the potential v(s) is complicated enough to bar the hope 
of exactly solving Eq. (|89|l in a simple way. Reasonable 
approximate solutions can be obtained by means of the 
variational method of Rayleigh and Ritz. We thus look 
for the minimum of the average energy 



(it, H u) 
(u, u) 

whose reduced and dimensionless form is 

J Q oc [—w(s)w"(s)s — w(s)w'(s) + v(s)w 2 (s)s]ds 
f oc w 2 (s)sds 



(90) 



(91) 



The minimization can be carried out for different values 
of the coupling parameter 7, but the frozen core model 
is realistic only for small values of 7: a strong coupling 
with a frozen core would have the paradoxical effect of 
totally expelling the distribution function of the test par- 
ticle from the center of the beam, while keeping the Gaus- 
sian core always concentrated in the center. To elaborate 
an example we have chosen 7 = 3 and a test function of 
the form 



wis) = ^— (1 + as 2 + bs 4 ) 



(92) 



The values of the parameters that minimize the energy 
functional turn out to be 



a = 1.24; a = 0.63: 



0.03: 



2.92. 



(93) 
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FIG. 13: The variational approximation to the density dis- 
tribution produced by the space charge of a frozen core with 
7 = 3 (solid line). It is compared with the Gaussian solution 
(dashed line) corresponding to the case of a purely harmonic 
potential (7 = 0). 
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FIG. 14: Zoom up of Fig. ll3l that shows a node in the position 
density distribution at a » 4.7. 



metric beam we introduce the normalized radial density 
distribution 



p(r) 



A 



p 2 a 2 T{q + l) \2p 2 a 



(94) 



which is composed of a Gaussian core with variance a 
(the simple harmonic oscillator ground state) , plus a ring- 
like distribution whose size is fixed through the two pa- 
rameters p > and q > 0. The parameter < A < 1 
is the relative weight of the two parts. This density dis- 
tribution has the required form for suitable values of the 
parameters, but, at variance with the two previous exam- 
ples of approximate distributions (see Figs. 1111 and I14fl . 
it has no nodes. This is convenient for two main rea- 
sons: first because it is a rather general requirement for 
a ground state to have no nodes (this fact is a rigorous 
theorem for one-dimensional systems) . Moreover, it has 
been shown in Refs. [3(], H2, H3 that stationary distri- 
butions without nodes are also attractors for every other 
possible (non extremal) initial distribution: a property 
that will be useful in a future discussion of the possi- 
ble relaxation of the system toward a stable beam halo. 
For this cylindrically symmetric case the expressions <|30|) 
and H37fl for the velocity field and the potential become 



u(r) = yfpjr) 



; (+) 



(r) 



a u'(r) 
m u(r) 

„ 2 



V(r) 



a 

2m 



u"{r) lu'(r) 



2ma 2 2m \ u(r) r u(r) 
Moving to dimensionlcss quantities we have 



(95) 



The corresponding optimal radial density distribution 
sw 2 (s) is shown in Fig. ll3l where it is also compared with 
the Gaussian solution (dashed line) corresponding to the 
case of a purely harmonic potential (7 = 0). Notice that 
this density distribution has a form similar to that of the 
self-consistent solutions elaborated in Section ITVI for val- 
ues of 7 of about 10 (see Fig.|2J). It has a node which is 
not visible in Fig. ED We zoom up its plot in Fig. [21 
This behavior suggests that the halo distribution could 
be described, in first approximation, as a Gaussian core 
distribution plus a small ring of particles surrounding it 
and constituting the halo. Thus, in the following we will 
simply assume a specific form of a beam with a halo with- 
out deriving it as an effect of space charge interactions. 
Since it is not clearly established that a halo can be due 
only to space-charge effects, starting with a realistic ring 
distribution and trying to understand the dynamics that 
can produce it could be very useful in this respect. The 
techniques introduced in Section [IIII will be instrumental 
to this end. For a three-dimensional, cylindrically sym- 



w 2 (s) 



b(s) 
v(s) 



aV2' 



w 



2Ae 



2a 2 p{r) 



2(1 -A) 



p 2 T(q + l) 



crV2 



ma\/2 



a 

4ma 2 



v (+)( r ) 



a* 



aV2 

w"{s) + 1 w'(s) 
w(s) s w(s) 



V(r) 



p 

w'(s) 
w(s) 



(96) 



A three-dimensional plot of the dimensionless density dis- 
tribution is given in Fig. 1151 where we have chosen p = 1, 
q = 24 and A f» 0.49. In this example the importance 
of the halo ring has been exaggerated with respect to 
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FIG. 15: Plot of the cyli ndrically symmetric density distri- 
bution of the beam w 2 (v x 2 + y 2 ) given by Eq. I|9rip . which 
shows a halo ring surrounding the beam core. 




FIG. 16: The radial, dimensionless, forward velocity field for 
the distribution of Fig. 1151 (solid line), and for the ground 
state of a harmonic oscillator (dashed line) . 



the real case in order to make the effect clearly visi- 
ble. The explicit analytic expression of b(s) and v(s) is 
rather lengthy and not particularly illuminating: we sim- 
ply show the behavior of these functions, for the same 
values of the parameter for the previous example, re- 
spectively in Fig. and Fig. El where they are also 
compared with the corresponding dimensionless quanti- 
ties (— s and s 2 ) for the ground state of the harmonic 
oscillator. The function v(s) is the potential that in the 
stochastic model keeps the beam in the stationary state 
with a halo ring. 

In order to illustrate the main effects involved, let us 
give here the simpler formulae relative to the Id case. 
From now on we will consider only Id processes denot- 
ing by x one of the transverse space coordinates. We 
assume that the longitudinal and the transverse beam 
dynamics can be deemed independent, with the further 




FIG. 17: The radial, dimensionless, potential for the distri- 
bution of Figure ITS1 f solid line), and for a harmonic oscillator 
(dashed line). 



simplification of considering decoupled evolutions along 
the transverse directions x and y. Under these conditions 
the density distribution with a halo ring now reads 



p{x) 



A- 



-x 2 /2a 2 



OV27T 

+(1-A). 



( x 2 



so that the corresponding velocities and potential are 

u(x) 
v {+) (x) 



(97) 



a u'(x) 



777, u(x) 



V(x) 



a 



a 2 u"{x) 
Ama 2 2m u(x) 



+ 



(98) 



In terms of dimensionless quantities we then have for the 
distribution 



w 2 (s) 



aV2 



1 



2a 2 p(x) 
A) 



P 



pT (q 

while the forward velocity and the potential are 

x 



-o (99) 



b(s) 



v(s) = v 



aV2 



mu\f2 , . w'(s) 



a 

Ama 2 



1 



w(s) 

w"(s) 
w(s) 



(100) 



(101) 



In Figrs.^OEB and 1201 we respectively show the density 
distribution, the velocity field, and the potential for p = 
1,(7 = 24 and A = 0.85. As in the three-dimensional case, 
the importance of the halo ring has been exaggerated 
with respect to the real case in order to make the effect 
clearly visible. 
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FIG. 18: Plot of the Id density distribution @3 with a halo 
ring surrounding the beam core. 




FIG. 19: The dimensionless velocity for the Id distribution of 
Figure ITSl (solid line), and for a harmonic oscillator (dashed 
line). 




FIG. 20: The dimensionless potential for the Id distribution 
of Figure^] (solid line), and for a harmonic oscillator (dashed 
line) . 



VI. CONCLUSIONS 

In this paper we have presented a dynamical, stochas- 
tic approach to the description of the beam transverse 
distribution in the particle accelerators. In the first part 
we have described the collective beam dynamics in terms 
of time-reversal invariant diffusion processes (Nelson pro- 
cesses) which are obtained by a stochastic extension of 
the least action principle of classical mechanics. The dif- 
fusion coefficient of the process is identified with a lower 
bound for the omittance of the beam as discussed in the 
Section Hvl The collective dynamics of beams is then de- 
scribed by two nonlincarly coupled hydrodynamic equa- 
tions. This set of equations shows some formal analogies 
with the effective Ginzburg-Landau and Gross-Pitaevskii 
descriptions of the collective dynamics of self-interacting 
quantum many body systems in terms of classical, non- 
linear Schrodinger equation. We have shown that, in 
the framework of SM, it is possible to have transverse 
distribution which show a broadening and an emittance 
growth, typical of the halo formation. The interest of this 
approach lies in the fact that its dynamical equations al- 
low us, at least in principle, to determine the possible 
control potentials that are responsible for the stationar- 
ity of the halo. A few introductory examples of these 
potentials and distributions have been discussed in the 
section The controlling potentials can be engineered 
by suitable tuning of the external electromagnetic fields. 
In a previous paper |2Cj we have considered evolutions 
that drive the beam from a less collimated to a better 
collimated state, and we have furthermore shown that 
this goal can also be achieved without increasing the fre- 
quency of the betatron oscillations which can in fact be 
independently controlled during the evolution. In forth- 
coming papers we plan to elaborate on the extension of 
these techniques to the problem of engineering suitable 
time-dependent potentials for the control and the elimi- 
nation of the beam halo. In particular, it will be shown 
that the transition functions of the Nelson processes can 
be exploited to control these evolutions. We also plan to 
extend the analysis of the present paper to different, non 
Gaussian beam transverse distributions which could bet- 
ter describe the existence of halo losses in terms of large 
ratios of the maximum displacement to the RMS size of 
the beam. 



APPENDIX A: COUPLED FIELD EQUATIONS 

Starting from the Maxwell equations for the space-charge 
field felt by our particle in the beam 



V-E sc (r,i) = 



PscO, t) 



VxE sc (r,i) = -d t B sc (r,t) 
V-B sc (r,<) = 

VxB sc (r,i) = (ioj ac (T,t) + j — - 
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and introducing the corresponding electro-magnetic po- 
tentials (A sc ,$ sc ) from 

B sc (r,i) = V x A sc (r,t) 

E ac (r,t) = -d t A sc (r,t)- V$ sc (r,t) 

with the gauge condition 

V-A sc (r,t) + la t $ sc (r,i) =0, (Al) 

cr 

through the usual procedure we get the wave equations 
1 



V 2 A sc {r,t)- -dfA sc (r,t) -- 
VH sc (r,t) - \d^ sc {r,t) = - 



MojscM) (A2) 
Psc(r,t) 



<~0 



(A3) 



It is also well known (see for example [T^, chapter XV) 
that the Schrodinger equation for spinless, charged par- 
ticles in an electro-magnetic field (A, $) is 



iad t ip 



1 

2m 



go* 



0. 



For a charged particle in the beam the electro-magnetic 
field is the superposition of the space-charge po- 
tential (A sc ,$ sc ) plus an external, control potential 
(A ext , <&ext), and hence we obtain the Schrodinger equa- 
tion 143|) that we rewrite here for convenience: 



1 

2m 



tQ; V- *(A S 
c 



A ea; {) 



(A4) 



It is apparent now that l|Aljl . (|A2jl . (|A3|t and (|A4jl con- 
stitute a self-consistent system of coupled, non linear dif- 
ferential equations for the fields ijj, A sc and $ sc . 

Since we are in a reference frame comoving with the 
beam, in the following we will suppose to deal only with 
stationary wave functions of the form (|42|) . In this case, 
as already remarked, we get j sc = 0, so that we can 
take the trivial solution A sc = of (|A2|) . The gauge 
condition 1|A1[) then implies that dt$ S c = and the wave 
equation (|A3|I reduces itself to the Poisson equation for 
the electrostatic potential. Finally, by supposing also 
A ex t = 0, our system is reduced to only two coupled, 
non linear equations for the couple (it, & S c), namely 



Eu 



Nqo \u\ 2 



eo ||w|| 



where we have defined 



\u\\ 2 = 



R 3 



\u(r,t)\' d A v 



Then if we pass to the potential energies by putting 
V ext {r, t) = q $ext{r, t) , V sc (r 7 1) = q <S> sc {r, t) 



our equations are reduced to 



— V 2 u +(E- V ext - V sc ) u = 0, 
2m 



Nq 2 \u\ 



e o \\u\\ 



(A5) 
(A6) 



Remark that here the wave functions ip and u, while cer- 
tainly normalizable, are not supposed in general to be al- 
ready normalized. However the equations (|A5|) and l|A6(l . 
albeit non linear, are apparently invariant for the mul- 
tiplication of u by an arbitrary constant. The equa- 
tions (IA5I) and (IA6I) will be used in the Section llVl 



APPENDIX B: CYLINDRICAL HARMONIC 
OSCILLATOR 

In the model discussed in this paper we suppose that 
the external potential V e is a cylindrically symmetric, 
harmonic potential with a proper frequency o>, so that 
the ground state without space-charge interaction will 
have a variance 

o a 



2mu! 



This means that, in the usual cylindrical coordinates 
{r, <p, z} with r = x 2 + y 2 , we have 



V e {r) = ^u 2 r 2 



and the corresponding Schrodinger equation 



(Bl) 



iadfip = -— V 2 V' + T4V' 
Zm 



a 
2m 



2 r 



d r (rd r ) + -d 2 +d 2 



ip + V e ip 



with periodic conditions at z = ±L/2 (for a bunch of 
length L) has the eigenvalues 



E nk = {n + l)auj + k 2 



2tto- 



and the normalized eigenfunctions 

i) nk l(r, y, z, t) = R nl {r) $,(¥>) Z k {z) e ~ iE ^ a 

with 



Rni(r) 
Z k (z) 



I mi 

Jltp 



aV2 



(I) 



2a 2 



2tt 

zi2k77z j L 



18 



where L^\x) arc the generalized Laguerre polynomials 
and 



n = 0,1,2,... 

0, 2, 4, . . . , n if n even 

1, 3, 5, . . . , n if n odd 

k = 0,±1,±2,... 



We suppose that, by neglecting the space-charge inter- 
action and in a comoving frame of reference, the system 
will be correctly described by the ground state 



•0000 (r) 



-r 2 /4a 2 



associated to the eigenvalue Eqq = ato. 



(B2) 
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